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r~| ■ Abstract 

I , In the (M+1)SSM an additional gauge singlet Weyl spinor appears in 

^j' the neutralino sector. For a large part of the parameter space this approx- 

,-S^ ■ imative eigenstate is the true LSP. Then most sparticle decays proceed via 

an additional cascade involving the NLSP -^ LSP transition, where the 

KJf I NLSP is the non-singlet next-to-lightest neutralino. We present a compre- 

^ ■ hensive list of all processes, which contribute to the NLSP — > LSP transi- 



a 



tion, the partial widths and the total NLSP decay rate. We perform a scan 
of the parameters of the model compatible with universal soft terms, and 
find that the NLSP life time can be quite large, leading to macroscopically 
displaced vertices. Our results imply that the signatures for sparticle pro- 
duction in the (M-l-l)SSM can be very different from the MSSM, and are 
important for calculations of the abundance of dark matter in this model. 



1 Introduction 

The supersymmetric extension of the standard model with an additional singlet 
superfield [^|-|^ has some attractive features: the superpotential can be chosen 
to be scale invariant, hence the only dimensionful parameters - and thus the 
electroweak scale - enter via the soft supersymmetry breaking terms. With a 
scale-invariant superpotential, and assuming universal soft terms at the GUT 
scale, the model has the same number of free parameters as the MSSM. Several 
analyses of the parameter space of the model have previously been performed 
in [§,§]. It has been found that a considerable region is consistent both with 
theoretical constraints (correct SU{2)l x f/(l)y symmetry breaking, no squark 
or slepton vev's, neutral LSP) and experimental lower bounds on sparticle and 
Higgs masses. 

It is very important to investigate in what respect the phenomenology of the 
(M+1)SSM differs from that of the MSSM. The signatures for sparticle produc- 
tion could be different, and one would like to know which processes could serve 
to distinguish the two models. 

The particle content of the (M+1)SSM differs from the MSSM in the form of 
additional gauge singlet states in the Higgs sector (1 neutral CP-even and 1 CP- 
odd state) and in the neutralino sector (a two component Weyl fermion). These 
states are mixed with the corresponding ones of the MSSM, and the physical 
states have to be obtained from the diagonalization of the mass matrices in the 
corresponding sectors. An interesting result of the analyses in [^, ^] is that, for 
most of the parameter space, the mixing angles involving the singlet states are 
actually quite small. Consequently there exist physical quasi singlet states which 
have only small couplings to the gauge bosons and the MSSM sparticles such as 
charginos, squarks and sleptons. These states then have only small production 
cross sections and it seems to be nearly impossible to observe them in present or 
future experiments. 

A notable exception can occur, however, in the neutralino sector. In the 
MSSM the neutralino sector consists of two gauginos (the bino and the neutral 
wino) and two higgsinos. Typically the LSP - the lightest supersymmetric par- 
ticle, which is stable if one assumes, as we do, R-parity conservation - is the 
lightest eigenstate of the neutralino mass matrix. The LSP will appear as one 
of the final states of each sparticle decay, and its non-observability is responsible 
for the well-known missing energy/momentum signature of sparticle production. 

The situation in the (M-l-l)SSM, where a singlino state is added to the neu- 
tralino sector, depends crucially on its mass with respect to the MSSM LSP mass: 
If the singlino is heavier, it will very rarely be produced and it will be practically 
unobservable. If it is lighter (and is thus the true LSP), it will now appear at the 
end of the decay chain of sparticles decays. To be more specific, from the analy- 
ses performed in ^ one finds that the MSSM LSP, within the allowed parameter 
space of the (M-l-l)SSM, is essentially a bino. In the singlino LSP case of the 
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(M+1)SSM one has to keep in mind the small couplings of the singlino to all the 
other particles. If the sparticles are heavier than the bino (which turns out to be 
always the case, except for some large supersymmetry breaking terms that yield 
sparticles out of reach for LEP2) they thus prefer to decay into the bino to which 
they couple more strongly. Only then the bino will decay into the singlino LSP, 
which will give rise to an additional cascade in the sparticle decay chain. Since 
this process modifies the signatures for sparticle production considerably, we will 
investigate the bino to singlino transition in detail in this paper. 

Production and decay of neutralinos have previously been discussed in the 
MSSM in, e.g., ref. and in the (M+1)SSM in ref. 0. (Many of the formulas 
of the partial decay widths in our appendix D can be found in these papers.) 
In ref. [^ production cross sections and branching ratios of neutralinos in the 
(M+1)SSM have been presented for several scenarios concerning the low energy 
parameters. 

Here, however, we are interested in a comprehensive analysis of the part of the 
parameter space of the (M+1)SSM which is compatible with universal soft terms 
at the GUT scale, which corresponds to the case of a singlino LSP, and where 
sparticle production is kinematically possible at LEP2. Our aim is to see, which 
bino lifetimes and bino branching ratios are possible under these assumptions. 
With our results at hand one can decide, which signatures for sparticle production 
(beyond the ones of the MSSM) are most promising in the framework of the 
(M+1)SSM, and which part of the parameter space can be tested. 

Our approach follows closely the one of |^, up to slightly different experimen- 
tal constraints on the parameter space: We start to scan the complete parameter 
space of the model (the universal soft terms and Yukawa couplings, see the next 
section) as defined at the GUT scale. For each point in the parameter space we 
compute the effective parameters at low energy by integrating the renormaliza- 
tion group equations from Mqut down to Mz- Then, for each set of parameters, 
we minimize numerically the effective potential including the one loop radiative 
corrections induced by top quark and squark loops. We check whether the ab- 
solute minimum of the potential breaks SU{2)l x f/(l)y as desired, whether 
squarks or sleptons do not assume vev's, and whether the physical masses of the 
top quark and the sparticles satisfy the (model independent) present experimen- 
tal constraints. (Details are given in the next section.) Finally we require that 
the LSP, Xi; is essentially a singlino state (otherwise the signatures for sparticle 
production are not different from the MSSM), and that the mass of the NLSP 
X2 (which is the bino) is below Mz- Under this approximate condition sparticle 
production at LEP2 is kinematically allowed. 

For each of the ~ 10^ remaining points in the parameter space, we compute 
the following decay widths: X2 ~^ Xi^^^~) X2 ~^ Xi^^y X2 ~^ XiQQ (i^i all cases 
we take Z, Higgs and slepton or squark exchange into account), X2 ~^ Xi+Higgs 
and X2 ~~* Xi + 7- (The radiative decays into a photon have previously been 
considered in the MSSM in, e.g., refs. ^,^, and in the (M+1)SSM in ref. 



[0.) The results give us both the hfe time and the branching ratios of the 
NLSP -^ LSP transition for each point in the parameter space of the model, 
which is consistent with universal soft terms, present experimental constraints, 
and which is of potential phenomenological relevance for LEP2. 

Clearly many steps of this procedure (e.g. the integration of the RGEs, the 
minimization of the effective potential, the diagonalization of the mass matri- 
ces and the phase space integrals) require numerical methods. These allow us, 
however, to obtain the results with satisfactory accuracy. On the other hand, 
we find it very desirable to understand at least the rough features of our results 
(and of the range of the low energy parameters) using analytic approximations 
to the integrated RGEs, the minimization of the effective potential, and the di- 
agonalization of the mass matrices. Therefore we spend some time in sect. 3 in 
order to discuss the interplay between the different theoretical and experimental 
constraints on the parameters within such analytic approximations. These ap- 
proximations allow us to understand the relative importance and the orders of 
magnitude of the different decay widths in sect. 4. The results on the different 
decay widths, branching ratios and the total life time presented in sect. 4 (and 
in the figures) are, however, based on the "exact" numerical procedure. 

Our results show that, even in the limit of tiny couplings of the singlino, a 
priori a large number of different processes can contribute to the bino to singlino 
transition. Only after a detailed investigation of all the partial widths we find 
that only a few of them are relevant: Essentially the three body decays with 
two leptons in the final state (via virtual slepton exchange) or with qq in the 
final state (via virtual Z exchange), and in some cases the two body decay into 
a real singlet Higgs scalar or a photon. Interestingly enough we find that, for 
small enough Yukawa couplings, the lifetime of the bino becomes so large that 
displaced vertices appear to be visible. 

Two cosmological issues should also be discussed within the (M+1)SSM, 
namely domain walls and dark matter. The (M+1)SSM with a scale invariant 
superpotential has a discrete Z^ symmetry, which can lead to the formation of 
domain walls with an unacceptable energy density during the electroweak phase 
transition 0. As discussed in p, possible ways out of this problem are to embed 
the discrete symmetry into a gauge symmetry at some large scale, or to add tiny 
mass terms, which do not modify the phenomenology in a visible way, but which 
break the Z^ symmetry sufficiently such that the domain walls are removed. 

The LSP of any supersymmetric theory with conserved R parity is a priori 
a welcome candidate for cold dark matter. It will necessarily be produced in 
sparticle decays in the early universe, and its relic density will strongly depend 
on its annihilation cross section. The (M-l-l)SSM has been considered in this 
respect in |Q], where both upper and lower limits on the LSP relic density have 
been imposed. In [Q it has been argued that the singlino LSP scenario of the 
(M+1)SSM is essentially ruled out, since the pair annihilation cross section is too 
small, and consequently the relic density is too large. However, in [Rl only the LSP 



pair annihilation has been considered. In particular in the case of small Yukawa 
couplings the situation for a singlino LSP is, however, much more complicated: 
The binos could pair annihilate before the LSP is produced, and the bino-singlino 
coannihilation rate is much larger than the singlino pair annihilation rate. In 
order to determine the dark matter constraints in the (M+1)SSM with singlino 
LSP reliably it is thus absolutely necessary to know the bino lifetime or the bino to 
singlino decay rate. Apart from the modified signatures for sparticle production 
the results of this paper will thus also find applications in the investigation of the 
dark matter in the (M+1)SSM. 

The paper is organized as follows: In the next section we present the la- 
grangian and discuss briefly the method of the scanning of the acceptable param- 
eter space; this procedure follows the one of ref. p. In section three we study 
the range of parameters in the singlino LSP scenario in some detail, focussing 
on analytic approximations. In section four we investigate all possible bino to 
singlino decay channels, the corresponding contributions to the partial bino decay 
widths, and the bino lifetime. We present both approximate analytic results, and 
"exact" results based on the numerical procedure. In section five we discuss our 
results and its phenomenological consequences. 

2 Parameter Space of the (M+1)SSM 

In this section, we study the parameter space of the model with the same as- 
sumptions as in p. The superpotential of the (M+1)SSM is given by 

W = \H1H2S + \kS^ + ... (1) 

where the ellipsis stand for quarks and leptons Yukawa couplings, 

Hi={%\. H2=i^^\ and H^H^ = H^H, - H^.H^ (2) 

Here the Higgs doublet Hi couples to the up-type quarks, and H2 to the 
down-type quarks and the charged leptons. Therefore the usual parameter j3 is 
given by 

tan/? = — (3) 

with hi = (H^). S denotes the gauge singlet superfield beyond the MSSM. The 
superpotential contains no HH1H2 term. An effective /i term is generated once 
the scalar component of the singlet 5* acquires a vev s: 

/i = As. (4) 



The only dimensionful parameters of the model are the supersymmetry break- 
ing gaugino masses, scalar masses and trilinear couplings (for simplicity we do 
not display the terms involving squarks or sleptons): 

Csoft = ^(M3A^A^ + M2A*2A^2 + ^iAiAi)+ h.c. 
-ml\Hi\^ - ml\H2\^ - ml\S\^ 
-\AxH1H2S - \kA^S^ + h.c. (5) 

Ai, A2 and A3 are the gauginos of the U{1)y, SU{2)l and 5'f/(3)c gauge groups 
respectively. The scalar components of the Higgs in @ are denoted by the same 
letters as the corresponding chiral superfields. 

The scalar potential contains the standard F and D terms, the supersymmetry 
breaking terms and one loop radiative corrections of the form 

1 { M'^\ 

t^-=64^ST.A<^>n(^). (6) 

In (^) we take into account only top quark and squark loops, but we include 
the numerically important contributions beyond the leading log approximation 
which result from the complete top squark mass matrix. Q^ denotes the renor- 
malization point, and all the parameters in eqs. (|I]), (^ and (^ have to be taken 
at the scale Q"^ ~ M"^. 

The supersymmetry breaking terms of the model are constrained by requiring 
universal terms at the scale Mqut ~ 10^^ GeV. The independent parameters of 
the model are thus universal gaugino masses Mq (always positive in our conven- 
tion), a universal mass for the scalars ttiq, a universal trilinear coupling Aq (either 
positive or negative), and the Yukawa couplings Aq and kq of the superpotential 
(|ID at the scale Mqut- In addition the top quark Yukawa coupling affects the 
renormalization group evolution of the parameters from Mqut down to the elec- 
troweak scale. The value of the Z mass fixes one of these parameters with respect 
to the others, so that we end up with 5 free parameters at the GUT scale, as 
many as in the MSSM with universal soft terms. 

Following the same procedure as in 0, we perform a scan over the complete 
parameter space of the model at Mquti integrate the renormalization group 
equations (RGE) down to the electroweak scale, and minimize the low energy 
effective potential including the radiative corrections (^ numerically in each case. 
We check whether we have found the absolute minimum of the potential, and 
verify whether squarks or sleptons do not assume vev's, which would break color 
and/or electromagnetism. Already at this stage, the condition to avoid selectron 
vev's (which are the most dangerous ones) yields a constraint on the parameter 
space 0: 

^ ^ - 2.5 (7) 

Mo ^ ' 



In the remaining cases we diagonalize numerically the mass matrices, compute the 
physical masses of all particles and impose the following experimental constraints: 

m^, > M^/2GeV [12], 

mt = 175±6GeV ^. (8) 

Note that, since signatures for sparticles production in the present scenario 
may be different from the MSSM, we cannot apply the standard MSSM analysis 
to the latest data from LEP1.5 and LEP2. This data should rather be reanalysed, 
in the context of the (M+1)SSM, using the results of the present paper. However, 
the LEPl results on the Z width and thus the sneutrino mass mc, remain valid. 
Moreover it turns out that the essential properties of the neutrahno sector do not 
depend on the details of the lower limits on, e.g., the chargino or slepton masses. 

Furthermore, within the present assumption of universality of the soft pa- 
rameters at the GUT scale and the singlino LSP scenario, eqs. (|]) imply already 
strong constraints on the other new particle masses (cf section |^), so that nearly 
all the other experimental bounds turn out to be automatically satisfied: 

m^o+m^o>Mz or r(Z -^ x-X?) < 7MeV if {t,j) = (1,1), 



r(Z^X-X-)<30MeV if (2,j)7^(l,l), 



mH± > 150 GeV , rriAo > 130 GeV, 
rrii^ > 190 GeV , m-g > 280 GeV, 
m + > 60 GeV , m,- > 60 GeV 

where y4° is the lightest non-singlet neutral CP-odd Higgs. The lightest non- 
singlet CP-even Higgs is in the range from 100 GeV up to 140 GeV. 

As emphasized in |Q,|§], the allowed parameter space of the (M+1)SSM is in 
general characterized by small values of the Yukawa couplings A and k 
(A, K < 0.3). As we will see in the next section the singlino LSP case corre- 
sponds to even smaller values of the Yukawa couplings A, k < 10~^. 



3 Singlino LSP Scenario 

In this section we present some special features of the singlino LSP scenario. 
In particular, we derive some approximative constraints on the high energy free 
parameters and some analytic approximations for the low energy masses and 
mixing factors. These approximations are useful to understand the features of our 
numerical results and will provide us with helpful guidelines for the calculations 
of the next section. To this end we use approximate relations between the low 
energy and high energy parameters of the model (as obtained from the RGEs), 
and relations obtained from the minimization of the tree level potential. At the 



beginning we assume that the Yukawa couphngs A and k are small, but this 
assumption will be justified below. 

In order to derive the constraints on the parameters implied by the singlino 
LSP scenario we first have to find approximate expressions for the singlino and 
the lightest non-singlet neutralino masses. 

Let us start with the singlino mass M. From the neutralino mass matrix 
( |6lD of Appendix ^ one finds that the mixing of the singlino to the higgsinos 
is proportional to A and thus relatively small. Hence, the singlino remains an 
almost pure state, and its mass is 

M = 2ks. (9) 

Using the minimization of the tree level potential the vev s and hence M can be 
related to the bare parameters of the model: For small Yukawa couplings (and 
hence s ^ /ii, /i2) the minimization equation (|69D for the singlet becomes 





(10) 

A^ and ms being only slightly renormalized between Mqut and Mz (cf Ap- 
pendix P), one obtains the singlino mass in terms of the GUT parameters: 

(11) 

Note that M has the sign opposite to Aq. The condition for the minimum (^) 
to be deeper than the symmetric one {hi = h2 = s = 0) reads 

Al > 9ml (12) 

so that 

^l^ol < \M\ < \Ao\. (13) 

Next, we estimate the masses of the lightest non-singlet neutralinos. The 
higgsino effective mass term /x = As turns out to be quite large (see below). 
Since the mixing terms between the gauginos and the higgsinos are at most of 
0{Mz/2fi), one finds that, to a good approximation, the hghtest non-singlet 
neutralinos are the (nearly pure) bino and the wino. Their masses Mi and M2 
are related to Mq as given in appendix j^: 

Ml ~ .5M2 ~ .41Mo. (14) 

The condition for the singlino to be the LSP is given by M < Mi, which, 
combined with (|T3|) and ([T^), yields 

1^1 < .6M0. (15) 



(Note that this condition is compatible with the necessary condition for the ab- 
sence of color and/or electromagnetism breaking vevs, eq. (0).) Eq. ([TsD together 



with eq. (12) implies that the singlino LSP scenario discards large |y4o|/Mo and 
itiI/Mq ratios, and it just corresponds to a very natural "gaugino dominated 
scenario" , gaugino masses being the largest soft terms. Then, the masses of all 
non-singlet sparticles can be expressed in terms of Mq and are therefore strongly 
correlated. 

For later use it is convenient to introduce a parameter rj, defined by the ratio 
of the masses of the lightest and next-to-lightest neutralinos: 

m^o M 
m^o Ml 

Unlike in the MSSM, it is not fixed by universality constraints at the GUT scale, 
but it is rather a free parameter varying from —1 to +1. Eqs. ([T3| ) and (|14[) allow 
us to express t] easily in terms of the bare parameters Aq and Mq: 

Next, we wish to estimate the higgsino effective mass term fi = Xs, and show 
that it is quite large. Below, the knowledge of /i will allow us to relate r] to 
the Yukawa coupling A. First, the minimization of the tree level potential with 
respect to hi and /12 gives the relation 

,an';3.4±4±i£|Z?. (18) 

mf + /i2 + M|/2 ^ ^ 

The condition for a non-trivial minimum corresponds to the condition that the 
denominator of eq. (|l^) has to be positive. The approximative solutions of the 
RGEs ( P^ and ( P3| ) imply that m| > 0, whereas ml < 0. Hence one obtains 



f^"> -mi- .5Mi > 2.IM0' - .5M^. (19) 

Since, in addition, phenomenological constraints imply that Mq is quite large 
(Md ^ 110 GeV, see below), one finds /i^ 3> M|. Actually, from our numerical 
results (within the singlino LSP scenario) we find tan/? > 6 which implies 

fx^ ~ 2.5M^ - .5M|. (20) 



in agreement with (|T^). 

Next, we derive an upper limit on the Yukawa couplings within the singlino 
LSP scenario. First, from the absence of a deeper unphysical minimum of the 
Higgs potential with h2 = s = 0, the following inequality can be derived P]: 



«:< 3.10-4- (21) 



-2 ^0 
Mi 



Furthermore, from the numerical analysis, the Yukawa couplings A and k turn 
out to be closely related: 

K ~ A^-^±-l (22) 



Using ([TTD in (^) and (|2^) one finds that the singlino LSP scenario requires small 
Yukawa couplings, X, k, < 10~^. 

On the other hand, no lower limit on the Yukawa couplings has been found 
in our analysis; we allowed for couplings as small as A = 10^^. In this regime one 
can show that the singlino LSP scenario follows automatically: Using @, (^, 
and pop, one gets 



\M\ = 2k\s\ = 2fi;|/i|/A ~ 2A-^^-V| ~ 3.2A-^^-^Mo, (23) 

so that 

\rj\ ~ 7.7A-^±-l (24) 

Hence, for very small values of A (A < 10~^ in our numerical analysis), the singlino 



LSP scenario is always realized and r^ ^ 1. The compatibility of eq. (24) with 
eq. (p!7| ) requires some relation between the bare parameters Aq, Mq and Aq of 
the model which is, however, not very stringent. 

Herewith we conclude the discussion of the constraints on the parameters of 
the model implied by the singlino LSP scenario. With the help of these results we 
can now obtain approximate expressions for all quantities which are required in 
order to calculate the bino to singlino decay widths: the mixing parameters of the 
singlino with the other non-singlet neutralinos, and the masses of the sfermions 
and the Higgs bosons. Let us first study the mixing parameters: 

For \r]\ not too close to 1, one can expand the singlino and the bino eigenstates 
in small mixing parameters in the basis of eq. (E9|) : 



N- 



ij 



Xg2{hl - hl)fi Xgi{hj - hl)fi 



\V2{M - M2)(M2 - /i2) ' y2(M - Mi)(M2 - /i2) 
A(/i/ii - M/ia) A(/i/i2 - Mhi) \ 



2 



M^-fi^ ' M2 - /i' 
( -gig2{Mi{hl + hi) + 2fxhih2) gijiih^ + M^hi) 



(25) 



2(M2-Mi)(M2-/i2) ' ' y2(M2-/x2) 
-gi{Mih2 + fihi) Xgiihf - hl)fi \ 



v^(M2 - /i2) ' v^(Mi - M) (M2 - /i2 



(26) 



Using eqs. (piil) , (pO|), ([16|) and /ii ^ /12, these components can be expressed in 
terms of 77, Mi and Mz- (However, (p5D and (^) are not valid anymore in the 
degenerate case |?7| -^ 1.) These expressions for the mixing parameters will be 
used extensively for the analytic approximations in the next section. 
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Next, we turn to the sferinion sector. The hghtest states are the sneutrinos 
v and the "right handed" charged sleptons Ir. The approximate expressions for 
their masses are (cf Appendix [C]) 

m?^ = m| - sin^ OwM^ cos 2(3 ~ IhMl + .23M|, (27) 

2 - -2 ^ iM| cos 2/5 ~ .b2Ml - .5Ml (28) 



rrij^ = m 



The lower hmit on rn-p (|^) combined with eq. (^) gives a lower limit on Mq 
( > 110 GeV), which in turn puts a lower limit on m^ ( > 60 GeV). "Left 
handed" charged sleptons and squarks are always much heavier, hence uninter- 
esting for the present phenomenology. 

The lower limit on Mq also yields a lower limit on the bino mass: 
m^o > 30 GeV. Subsequently we restrict our analysis to the regime m^o < Mz, 
where sparticle production at LEP2 (at least the pair production of binos) is 
kinematically allowed. In terms of Mq this corresponds to Mq < 230 GeV. In 
this region of the parameter space, the bino is the NLSP. Note that for larger 



values of Mq the D-term in eq. (p7|) becomes negligible and one has 



mf ~ .15M^ < Mf ~ .ITMq. (29) 



Thus, for Mq > 320 GeV, Ir turns out to be the NLSP. We shall come back to 
the case m^o > Mz in the last section. 

The large value of /i implies that the lightest chargino xt is mainly a wino of 
mass M2, which is related to Mq by eq. (p^. However, for small values of Mq the 
higgsino component can be quite large (up to 50%) and m ± smaller than M2. 
The lower bound on Mq yields m + > 60 GeV. 

Similarly, from the lower bound on Mq, we obtain the lower bound on the 
gluino mass given in the previous section. 

Finally, we briefly focus our attention on the Higgs sector. We see from the 
mass matrices given in Appendix |B| that the mixings of the CP-even and CP-odd 
singlets with the non-singlet Higgs fields are proportional to A, hence small. Here 
again, the singlet sector decouples and we end up with two almost pure singlet 
states, a scalar of mass 



-') 



^1,33 ^ tV^o - 8mg ( \Aq\ + ^Al - Smg ) , (30) 



and a pseudoscalar of mass 



3, 



^1,22 ^^l^oK 1^1 + \/^g-8m2) (31) 



where we have used eqs. ([72D, (fT^) and (p^O]). The pseudoscalar is always heavier 



than the scalar singlet and, using arguments similar to (^), one finds that their 
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masses are both roughly proportional to M ~ A'^^'^Mq in the singlino LSP case. 
Therefore, the singlet states are the lightest Higgs states. 

In the non-singlet sector, we have one CP-odd pseudoscalar A'^ of mass 

_2 0.3M2 ^ ^ 

where we have used (^) and the approximate solution of the RGE for Ax (PT|) 
in the limit mo,^o ^ Mq. As tan/5 is always quite large, A° turns out to be 
relatively heavy as already remarked in section |^. The mixing term between the 
CP-even fields iJi/j and H2R is proportional to: 

M^ -M^ -cot/3«l. (33) 

ifiij and if2ij are then almost pure states of mass 

M|ii ^ ^i (34) 

MI22 ^ 0.3Mo^tan/5>M|ii. (35) 

These approximations must be taken with care, as they do not include the nu- 
merically important radiative contributions beyond the RGE to the effective po- 
tential d^). Nevertheless, we find from the numerical analysis (with the radiative 
corrections to the effective potential included) the following particle assignments 
and mass ranges in the Higgs sector, in agreement with our rough estimates: 



^1 


r\^ 


Sr 





< 

r\-/ 


ms. 


< 

r\-/ 


60 GeV, 


^2 


r\^ 


HlR 


100 


< 

r\-/ 


ms. 


< 

r\-/ 


120 GeV, 


^3 


r\^ 


H2R 


130 


< 

r\-/ 


ms. 


< 

r\-/ 


360 GeV, 


Pi 


r\^ 


Si 





< 


mp^ 


< 


130 GeV, 


P2 


rvj 


A° 


130 


< 


nip^ 


< 


350 GeV. 



However, the upper bounds given above increase if one relaxes the condition 
Ml < Mz, allowing for higher values of Mq. 

4 Bino to Singlino Transition 

In this section we compute the bino to singlino decay widths. As already men- 
tioned, this process is crucial as it will appear at the end of every sparticle decay 
chain in the singlino LSP case. The different contributions are shown in figs. |l]-^. 
Exact formulae for the corresponding decay widths are given in Appendix]^. The 
production and decay of neutralinos have already been studied in the (M+1)SSM 
framework for a few selected points in the parameter space [^ . Instead, we have 
computed the partial and total decay widths numerically for each point in the 
parameter space obtained from our numerical scanning. 
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In the following, we first present some simple analytic approximations so as 
to understand the main features of the bino to singlino transition. Then we 
discuss our the "exact" results, which are based on the numerical integration of 
the RGEs, the numerical minimization of the full Higgs potential, the numerical 
computation of the mixings and mass eigenvalues in the neutralino, chargino and 
Higgs boson sector, and the integration of the exact phase space integrals. These 
results turn out to be in good agreement with the analytic approximations. 

First of all, let us consider the tree level three body decay X2 ~~^ Xiff oi fig. |l[ 
The fermions can be charged leptons, neutrinos or quarks (in which case we end 
up with two jets). All the decay widths are proportional to A^ - one factor A from 
the non-singlet component of the singlino, raised to the square - and are hence 
equally suppressed. Therefore, for each final state, we have to check whether the 
virtual Z, sfermion or Higgs exchange gives the main contribution and whether 
we have to compute interference terms. 

Let us start with a pair of charged leptons in the final state. The partial 
width via virtual Z exchange is given by eq. (^7^ [0|. It depends on the mixing 
factor O12 in the coupling Zx1x2 defined by eq. (0), and a phase space integral 
Iz defined by eq. (^8|). In our analytic approximations we assume a very light 
singlino, i.e. rj small. (We shall come back later to the case 77 7^ 0.) Using (p5|) 
and (|26| ) in the limit of large tan/3 and \rj\ -C 1, the mixing factor O12 can be 
written 

Oi2~1.710-2A -/ . (36) 



As we take \ri\ small, the phase space integral Iz is of 0(10^^), so that the decay 
width reads 

r(X2 ^ xTn - Q.10-'X'M^Iz{v,uJz) ~ lO-i^A^Mi. (37) 



For the slepton exchange, the partial width is given by eq. ( p.08| ) . Since the 



"left" type charged sleptons are always much heavier than the "right" type ones, 
their contribution will be relatively unimportant. The vertex factor involves the 
mixing factor N12 defined in eq. (p^). Using (|25| ) with the same assumptions as 



above yields 

\/2n.M.n ~ .13A ( 

Ml 



2 
V2(7iiVi2~.13A(^) . (3^ 



Mz 



The partial width then reads 

r(x° ^ xTn ^ 2.10-«A^Mi I ^ 1 Ij^iv,^lJ- (39) 



rriT 



As for the Z exchange, the phase space integral /;- , given by eq. (|109|) , is of 
0(10^^). One can infer from ( pTf ) that the ratio Mzjvni is always of 0(1). 
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Eq. (|39[) then gives 



Tixl^xTn ~ 10-^A^Mi (40) 

In the case of virtual Higgs exchange, the partial width is given by eq. ( |103| ) 
and depends, in this case, on the mixing factors Qau and Qai defined in eqs. (|99D , 
(|101|) , respectively. (Here and below the index / denotes a charged lepton and 
replaces the index / in eq. ( |101| ).) First, we observe that if the lightest Higgs 
scalar (which is the singlet) is too heavy to be produced on shell, the partial width 
for its virtual exchange is proportional to A^, and hence completely negligible: 

Qn2~A2 and Qn ~ A =^ t{xI ^ Xil^H - >^'- (41) 

The result is similar for the singlet pseudoscalar Sj, which is always heavier than 
Sr. As shown in the previous section, the second scalar 5*2 is mainly Hm, so that 

Q2i.-4iV24-|iV,3-^^, (42) 

miS22 mi 

'^'' " Tar " vf^ *''* 



where m; denotes the lepton mass and we have used eqs. (|25|), (p^) and (pSf). 
Eq. (|103|) then gives 



r(x° ^ xrn - io-^A^^/2 



niTT 



{V,^2). (44) 



As before, the phase space integral I2, given by eq. (|104|) , is of 0(10 ^). The 



only leptonic final state with sizable couplings to the Higgses is the t^t pair. 
Taking for Hifi a mass of order 100 GeV, we get 

ml ^ Xlr-r-) ^ 10--A^^^. (45) 



Even if one takes Mi = M^, this is completely negligible compared to (^01) . The 
second pseudoscalar A^ and the third scalar H2R being much heavier than Hm 
will give even smaller contributions. 

To summarize, the dominant contribution to the X2 ^ Xi^^^~ transition is the 
slepton exchange, and we do not need to compute any interference term between 
diagrams. This remains valid for any value of t], although the partial width can 
become significantly smaller than (|40| ) as the phase space is reduced for \ri\ -^ 1. 

Next, we turn to the neutrino production X2 ~^ Xa_^^- The Higgs exchange 
does not contribute. For the partial width via Z exchange, we get the same 
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result as for charged leptons with an extra factor 2 from the Zvv vertex. As the 
sneutrino is a "left" type sfermion, the vertex factor required for the sneutrino 
exchange is slightly different: 



2 



^(<72iVn-<7iiV,2)^.m(^) . (46) 



The partial width reads 



4 



> IQ-^X^Mi. (47) 

Although sneutrinos can be rather heavy (~ 150 GeV for Mi = Mz), this con- 
tribution always remains larger than the one from Z exchange. Thus, the virtual 
sneutrino exchange gives the main contribution to the X2 ~^ Xi^^ channel and 
the computation of interference terms is not needed. 

Finally we consider the decay into two jets X2 ~^ XiQQ- The top is too heavy 
to be produced. The partial width via virtual Z exchange is of the same order as 
in the case of leptons (with an extra color factor Ng = 3 and slightly different Zff 
couplings), whereas the squark exchange is strongly suppressed because squarks 
are always rather heavy. As for charged leptons, the virtual Higgs exchange plays 
no role. Hence, the virtual Z exchange is the only important contribution to the 
X2 ~^ Xa^QQ partial width, which therefore is always small compared to the partial 
width into two leptons via slepton exchange. 

With the approximate expressions for the three body decays X2 ~^ Xi^^^"; 
X2 ~^ Xi^^ aiid X2 ~^ XiQQ at hand, we turn now to the "exact" numerical 
results. In figs. ^^ we present our numerical results for the branching ratios of 
the three body decays X2 ^ X?^^^~, X2 "^ x\^^ and xl -^ x\m il = m, d, c, s, b) 
respectively, for ~ 10^ points in the parameter space described in sec. 2. Here we 
used exact expressions for the mixing factors, the phase space integrals and we 
included all contributions to a given final state. 

From the previous discussion, the branching ratios do not depend on A, since 
all the partial widths are proportional to A^, but essentially on the bino mass: 
For small values of m^o (~ 30 GeV), sneutrinos are lighter than charged sleptons 
(~ 45 GeV and ~ 60 GeV respectively), therefore the main contribution to the to- 
tal decay width is the neutrino production via virtual sneutrino exchange (fig. ^. 
As the bino mass increases, the sneutrino mass gets larger than the charged slep- 
ton mass. The dominant process is then the charged lepton production via virtual 
slepton exchange (fig. ^) (except for a small domain in the parameter space on 
which we shall come back in the next paragraph). As advertised earlier, the jet 
production (fig. ^) is always small. In the domain of large bino masses, where 
sleptons are also relatively heavy, the quark production via virtual Z exchange 
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can contribute up to ~ 20% to the total width. Genuinely we have 

Let us now study the two body decay into a real Higgs boson of fig. ||, start- 
ing again with approximate analytic expressions. The lightest non-singlet scalar 
(which is mainly Hm) is too heavy to be produced on shell. As already remarked 
in (0), one gets for the bino-singlino-singlet Higgs scalar vertex factor 

Qii2 ~ A^, (49) 

so that the partial width ( |11U| ) approximately reads 

Tixl^x'iSR)r^\'M,. (50) 

For small values of A, this is completely negligible compared to the three body 
decay rates. Yet, (|50|) involves no tiny numerical factor stemming from virtually 
exchanged particles. Hence, if A is not too small, real Higgs singlet production 
can even dominate the total decay width. However, as emphasized in the previous 
section, the masses of 5* and Sr are roughly proportional to A'^^'^Mq. Therefore, 
if A is too large, the singlet Higgs scalar and the singlino become too heavy to be 
produced on shell in the bino decay and this channel is kinematically forbidden. 

The numerical results are displayed in fig. |^. They are in good agreement with 
our approximations and one finds that for a small window in A, A ~ 10~^, the 
branching ratio of this process can reach 90%. In the same window, we could have 
the two body decay with a real singlet Higgs pseudoscalar X2 ~^ Xi^i- However, 
since the pseudoscalar singlet is always heavier than the scalar, this contribution 
remains small ( < 5%). If A > 2.10"^, the emitted Higgs singlet is heavy enough 
to decay into bb, which is then the main final state. For smaller values of A, 
this channel is kinematically closed. Depending on the singlet mass, the t~^t~ /cc 
channels are then favored. Smaller singlet masses correspond to smaller values 
of A, in which case the real Higgs singlet production is negligible. 

Finally, we turn to the radiative decay X2 ~^ Xil- A complete calculation 
involves loops with fermions + sfermions and charginos + W, charged Higgs 
and Goldstone bosons (depending on the gauge choice) |jTO|. The corresponding 
contributions decrease with increasing masses of the particles inside the loops. 
In the following analytic approximation, we then only consider the dominant 
diagram, involving the lightest particles in the loops, namely the "right" type 
charged sleptons (fig. ^). However, it is worth being stressed that we performed 
a complete numerical analysis, including all the loops mentioned above with the 
correct chargino and stop mass eigenstates |T^. The effective coupling ( |112| ) for 
three degenerate Iji loops is given by 
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Ml 



IliV.^n) 



(51) 



where /^, defined in eq. ( |113D , is of 0(10 ^) ii\rj\ -/^ 1. The partial width ( |111| ) 
then reads 



r(x^ - xh) - lo-^^A 



11X2^1 



Ml 



(52) 



Even for small values of Mi, this is totally negligible compared to the three body 
decay rates. This is not surprising since it is a contribution of higher order in 
perturbation theory. Note that there is no "dynamical enhancement" mecha- 
nism for this channel in our model as it can appear in the MSSM under special 
assumptions [10, 11|. However, there could be some "kinematical enhancement": 
Up to now, we have assumed \ri\ <^ 1 (i.e. very light singlino) in all our 
analytic approximations. What happens for |?7| ^ 1 ? On the one hand, all the 



three body decay phase space integrals (pS]), ( p.04| ) and (|109|) tend towards 0. As 
it has already been mentioned elsewhere [p!0|,p!l|, one can check that they are all 
of order 



i{v.^', 



H^i 



I'W 



(53) 



Hence, all the X2 ~^ Xiff channels are equally suppressed. Furthermore, since in 
this case the singlino mass is close to the bino mass, the two body decay with a 
real singlet Higgs boson is kinematically forbidden. 



On the other hand, it is well known that the radiative decay X2 



xh is 



usually less suppressed for 
around rj = ±1: 



\V\ 



1. One can expand the loop integral ( |113| ) 



I^{r],uJi 



ri^l 



Ti-^-1 



a + vf- 



(54) 



Therefore, the radiative decay gives the main contribution to the total decay 
width for rj -^ —1, but not for 77 — ;> 1. This phenomenon has not been observed 
before in the context of radiative neutralino decay. A similar effect exists for the 
neutrinoless double beta decay where the result depends on the relative sign of the 
Majorana neutrino masses |T9|. This rough estimate correctly fits our numerical 
results for the branching ratio Br{x2 —>■ Xi7); shown in fig. p. Actually, one finds 
that the main contributions to the radiative decay are the charged lepton/" right" 
slepton loops, the top/lighter stop loops and the lighter chargino/VT loops [ p^ . 
Interferences between chargino and sfermion loops being destructive, this leads 
to even smaller branching ratios for the radiative decay, and this channel is kine- 
matically enhanced only for a few points in the parameter space corresponding 
to ?7 ~ —1. 
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To conclude, we give in fig. ^ the total width of the bino to singlino transi- 
tion, including all the contributions discussed above, as a function of A. The A^ 
dependence is manifest for very small values of A where the singlino is always 
very light. As A increases, the singlino mass can take non negligible values, in 
which case the bino decay is kinematically suppressed. These cases correspond 
to the points in fig. ^ below the "fat" diagonal line. From the total width it is 
straightforward to compute the bino lifetime: 

h _ 6.58 10-25GeV.s 

For an energetic bino, the length of flight in the lab system is given by 

he 1.9710-i^GeV.m 



h% = V7' - 1 CT^l - ^.~o_,^oy^ - ^f^o^^o.^^ (56) 



One can then simply read off /^o from fig. |^. For T{x2 ~^ Xi^) ^ 10^^^ GeV 
(which corresponds to A < 5.10^^ or strong kinematical suppression), the bino 
escapes the detector, and the signature is the same as in the MSSM. In the 
other cases, we obtain the expected additional cascade, with a macroscopically 
displaced vertex (/^o > 1 mm) for T{x2 -^ Xi^) ^ 5.10"^'^ GeV. 

5 Conclusions and Outlook 

The purpose of the present paper was the calculation of the NLSP partial and 
total decay widths in the (M+1)SSM, in the case where the LSP is a singlino 
and sparticle production is kinematically allowed at LEP2. Then, the NLSP is 
essentially a bino and the bino to singlino transition appears at the end of all 
sparticle decay chains. (On the other hand, if the singlino is not the LSP, it 
will be nearly impossible to produce neither the singlino nor the Higgs singlet in 
collider experiments, since the singlet sector is always almost decoupled from the 
rest of the theory. The (M+1)SSM is then very difficult to disentangle from the 
MSSM.) 

We worked in the context of the constrained (M+1)SSM, with universal 
boundary conditions for the soft terms at the GUT scale. The essential features 
of this scenario have been discussed, using analytic approximations, in section 3. 
We have seen that, while the singlino LSP scenario is not a necessary consequence 
of the model, it corresponds to a natural "gaugino dominated" scenario, gaugino 
masses being the largest soft terms. Furthermore, the singlino is automatically 
the LSP for very small Yukawa couplings A and n. 

The gross features of the bino decay widths are easy to understand using 
analytic approximations. In section 4 we have presented and discussed such ap- 
proximate analytic expressions for the partial decay widths, which is in good 
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agreement with the results obtained numerically (without the corresponding ap- 
proximations). The numeric results have been presented in the form of the figs. ^ 
1^. (The corresponding detailed formulae are given in the appendices.) 

Let us summarize the behavior of the total decay width T{x2 -^ Xi^)- I^ 
principle, the bino partial and total decay widths can vary over many orders of 
magnitude, depending on the Yukawa coupling A, the bino mass m^o ~ Mi and 
the singlino to bino mass ratio rj (with rj ^ for A ^ 0). Generally one has 



r(x^^xW<10-^A^Mi. (57) 



For \ri\ <^ 1, the inequality (|57|) can roughly be replaced by an equality. For 
|?7| —^ 1, however, T{x2 — > Xi^) can be considerably smaller than the right hand 
side of ( p7|) because of kinematical suppression. The allowed range of the total 
decay width as a function of A is displayed in fig. |^, and the corresponding bino 
lifetime can lead to macroscopically displaced vertices. For very small values 
of r(x2 -^ Xi^): the bino may even decay outside the detector (in which case 
it imitates the true LSP of the MSSM). However, this requires tiny Yukawa 
couplings (A < 5.10"^) or strong kinematical suppression. Such scenarios could 
be probed by the same kind of apparatus as the slow neutralino to gravitino 
transition in the context of Gauge Mediated Supersymmetry Breaking models 



0- 

For a given value of T{x2 -^ Xi^)^ the branching ratios of the bino still vary 

essentially with the bino and singlino masses. In most of the parameter space, 

the three body decay (fig. |l]) dominates, and the relevant final states are uu, l^l~ 

or qq {q = u, d, c, s, b) and missing energy. For small values of the bino mass 

(~ 30 GeV), the uu channel dominates (fig §). Hence, the bino decays invisibly 

and its signature is just missing energy as for the true LSP of the MSSM. However, 

this channel never exceeds 90%, the remaining 10% corresponding to the visible 

l~^l~ channel (fig. ^. For larger values of the bino mass (up to Mz), on the 

other hand, the invisible final state ui? becomes less important and the charged 

lepton channel contributes up to 100%. The partial width into qq is always small 

compared to the partial width into l~^l~ , and we expect at most one jet event for 

ten charged lepton events. The characteristic signature for sparticle production 

would then be lepton events with high multiplicity (at least four, in e"*"e~ -^ X2X2) 

plus missing energy, eventually with displaced vertices. 

However, in the window 10^^ ^ A < 10^^, the two body decay X2 ~^ Xi'S'i 

dominates, if kinematically allowed (fig. |^. Si is then essentially the Higgs singlet 

with a mass varying between 3 and 35 GeV. If its mass is larger than ~ 10 GeV, 

5*1 decays into bb (with a branching ratio of ~ 90%); otherwise, t~^t~/cc are 

favored. (For such values of A, the bino and the real Higgs singlet will have short 

lifetimes.) The relevant final state would be two 6-jets (or eventually r+r^/c-jets) 

with an invariant mass peaked below 35 GeV. Such processes are totally exluded 

in the MSSM and would be a strong sign for the (M+1)SSM. 



Finally, in the degenerate case 77 ~ —1, all the previous tree level channels are 
kinematically suppressed, and the radiative decay X2 ~^ Xil dominates (fig. |]). In 
such a scenario, the bino would be very long lived (/^o > 1 m). (This corresponds, 
however, to a tiny fraction of the parameter space.) In contrast to the MSSM 



1^ , a dominance of the radiative decay is compatible with universal soft terms, 
and the (rather disfavored) condition tan/3 ~ 1 is not required in the (M+1)SSM. 

Herewith we have summarized our results, which have been obtained for 
m^o < Mz, the range accessible at LEP2. Let us, at the end, comment on 
the case of a bino which is heavier than the Z. One should consider two different 
regimes: 

In the intermediate range Mz < m^o < 130 GeV, the main decay mode 
becomes X2 "~* Xi^? with Xi and Z on shell. The total decay width is again 
proportional to A^, though larger than in the case of the three body decay (since 
no virtual particle needs to be exchanged). Hence, the bino would not be too 
long lived in this case, except for extremely small values of A. The characteristic 
signature for this additionnal cascade is missing energy plus the typical Z decay 
products. 

For a very heavy bino, m^o > 130 GeV, one has m^o > mj as already 
noticed in eq. (pOj). The "right" charged sleptons are hence light enough to be 
produced on shell, and the main channel is X2 ~^ ^^R- Since one needs at least 
m^o > 250 GeV in order to have m^o — mj > 10 GeV, the emitted lepton 
is very soft, hence difficult to detect. Yet, in this case, the true NLSP is the 
charged slepton and the process of interest is Ir — > /^Xi? appearing at the end of 
all sparticle decay chains. Then one obtains a charged slepton (eventually long 
lived, depending on A) decaying into an energetic lepton plus missing energy (the 
singlino). This case corresponds, however, to a very heavy sparticle spectrum, 
disfavored by solutions to the hierarchy problem. 
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Appendices 



A Neutralino Sector 



The mass terms for the neutrahnos are given by the following part of the la- 
grangian: 



C 



Wi 



W2 



+XsiIjiiIj2 + Xhiilj2'ips + Xh2ipiips - I^SljJslps 



+^MiAiAi + -M2A^A^ + h.c. (58) 

where the two component spinors Ai, A2, ipi, ip2 and ips are the supersymmetric 
partners of the B, W^, Hf, H2 and S bosons respectively. We introduce the 5 
component neutralino vector ||16|| : 

Then the mass terms read 

where the (symmetric) neutralino mass matrix M° is given by 



{-iXl,-iXi,iJi,ilj2,'ips 
-(^0)^M°(^0) + h.c. 



(59) 
(60) 



/ 



M' 



Mo 



Ml 



-fi'2/ii g2h2 



\ 



\ 



gih 


-gih2 


n 


V2 




V2 


-Xh2 







-Xhi 
2ns ) 



(61) 



The physical mass eigenstates are obtained by diagonalizing M.^ with a uni- 
tary matrix A^: 



m^oSij 



N- N- M° 



(62) 



m^o being the mass eigenvalues in increasing order of the neutrahno states: 



Ki 



N,,ij' 



ijYj 



1 5. 



(63) 



We take N real and orthogonal. Then some of the mass eigenvalues may be 
negative. Finally, one obtains the proper 4 component neutralino eigenstates by 
defining the Majorana spinors: 



X° 



Az 



1,...,5. 



(64) 
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B Higgs Sector 

We give here the potential, the minimization equations and the mass matrix 
for the neutral scalar Higgs without radiative corrections. The purpose of this 
appendix is only to set up our conventions and to provide some guidelines in order 
to make our analytic approximations easier to understand. A more complete 
analysis of this sector can be found in [Q. 

The scalar potential for the neutral Higgs fields is given by 

V = A2(|/70|2|5|2^|^0|2|^|2^|^0|2|^0|2)_^^2|^2|2 
-Xk{H',H^S*^ + h.C.) + ^(|i/iT - l^2T)' 



2 



_L^2| tt0\2 I ^2| tt0|2 , ^2 I q\ 
+m-^^\ll^\ +777,21X221 +771^|o| 

K A 

-XA^iH^HlS + h.C.) + —^{S^ + h.C.) (65) 



where g'^ = ^igf + g^)- We split the Higgs fields into real and imaginary parts: 

(66) 






The conditions for extrema of (|65f) are 



hi{\\hl + s^) + ^{h\ - hi) + m\) - Xh2s{Ax + Ks) = 0, (67) 

h2{\\hl + s^) - —{hi - hi) + ml) - Xhis{Ax + ^s) = 0, (68) 

s{\'^{hl + hi) + 2k'^s^ - 2\Khih2 + m|) - XAxhih2 + nA^s"^ = 0. (69) 



After the elimination of mf, ml and m| using eq. ( |67| - |69D , the elements of the 
3x3 symmetric mass matrix for the CP-even scalars in the basis (-ff?/j, -f^2_R5 ^r) 
are 

Mlu = g''hl + \s^{Ax + Ks), (70) 

til 

MI22 = g'hl + Xs^{A, + Ks), (71) 

"-2 

h h 
Ml^3 = \Ax^^ + Ks{A^ + 4:Ks), (72) 

^1,12 = i2X^-g^)hh2-XsiAx + Ks), (73) 

M|i3 = 2X^his + Xh2{Ax + 2Ks), (74) 

MI23 = 2X^h2S + Xhi{Ax + 2Ks). (75) 

Likewise, the elements of the 2x2 mass matrix for the CP-odd pseudoscalars 
in the basis {A^, Sj) where the would-be Goldstone boson has been projected out, 
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read 



Ml,,, = Xs{A, + Ks)^±^, (76) 

hih2 



t2 M a . a N^1^2 



M^22 = A(AA + 4fi:s)^^-3/tA^s, (77) 



M|i2 = A(AA-2/ts)y'/if + /ii. (78) 

The mass eigenstates of the scalars are denoted by Sa=i^2,zijnsi < ""^52 ^ ""^53) 
and those of the pseudoscalars by Pa=i,2ijnpi < '"^P2)- 

C Approximate Results of the Integrated RGEs 

In this appendix we display some simple analytic results of the integrated RGEs 
in the approximation where the dependence on all Yukawa couplings but the one 
of the top quark, ht, are neglected. Such solutions have been first discussed in 



the MSSM framework in ^\. The Renormalization Group Equations for the 
(M+1)SSM can be found in 0. We assume universality for the soft terms at 
the GUT scale and no flavor mixing. For a more complete and general set of 
solutions, cf iTH . 



First, let us define the parameter p by 

P = TT- (79) 

'''crit 

where hcnt is the infra-red fixed point solution for ht-, hcrit — 1-13. We find from 
our numerical results with a singlino LSP: -7 < p < .8. 

The Yukawa couplings A and k, are only slightly renormalized: 

A ^ Ao , K ~ kq. (80) 

The results for the soft trilinear terms read 

^A = Ao(l-0+(l.llp-.59)Mo, (81) 

A. = Ao. (82) 

Soft scalar masses are as follows 

m| = ml, (83) 

2 _ A 3 \ 2 P(l~P)/4 ooo/i#\2 



+ (.52 - 3.71p)M2, (84) 

ml = ml + .52M^, (85) 

ml = ml + .52M^, (86) 

ml = ml + .15Ml (87) 
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Finally, we have the usual gaugino mass relations 

M2 = .82Mo, Ml = ~M2 - Im2 = AlMo. (88) 

85-2 2 

D Decay Widths 

We first repeat some general features of the three body decay of X2 with mass m^o 
and momentum p2 into Xi with mass m^o and momentum pi and two massless 
fermions / and / with momenta k and k' respectively. 
Using the Mandelstam variables 

S = {p2-Pi)\ t = {p2-k)\ u={pi-k)\ (89) 

the differential decay width can be written as 

dV = f-^ y IMl^dudt (90) 

5127r3m3o ^^^ 

where A4 is the invariant amplitude for the processes under consideration and 
Nf is the color factor of the fermions. The different diagrams are shown in fig. |I[ 
As seen in the main part of the paper, we only need the integrated width for 
each process separately, without interference term (except in the case of sfermion 
exchange, see below). The integration limits are 

0< s < (|m^o| - \m^o\y 



'^min,max 



-2 H§ - "^x? - ^ ± .yA(mJo,mJo,5)) (91) 

u + t + s = mlo + mlo 

A2 Al 

A being the usual triangle function: 

A (a, 6, c) =a^ + b'^ + c^ - 2ab - 2ac - 2bc. (92) 

To begin with, let us consider the Z exchange. We use the following notations: 

• The parameters in the Z-fermions coupling are 

Lf = Tsf — Qf sin^ 6w, ^f = ~Qf ^in^ Ow (93) 

where T-^f and Qf are the isospin and the charge of the fermion respectively. 

• The Z propagator is given by 

Dz{x) = {x-Ml + tMzTz)-\ (94) 
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• The parameter in the ZX1X2 couphng is 

O12 = iVi3iV23 - NuN24 (95) 

where Nij denote the mixing components of the neutrahnos as given by 
eq. (H). 

The invariant amphtude of the process is [^ 

Y.\Mz\' = 4g\L} + R})0l,\Dz{s)\'x 

spins 

X ((^^0 —t){t — m^o) +t ^^ u + 2m^om^os] . (96) 

Hence, the partial width reads 

w~o z .Off, Nfg\L) + R))Ol, ml, 

r(X2 -^ Xjf) = f^ ^/z(r/,^z). (97) 



Iz is a phase space integral depending on r] defined in eq. (pig) and ujz = m^o/Mz'- 



1 f^-H A(42-l-2r7 + 3r/^)(l + r/-^) 
'^^^' ^^) = 3 h-f ^' {u\{2z^rf-V)-Xy ^'') 



with A = ^(1 - zf - if. 

Next we turn to the Higgs exchange. We only display formulae for a CP-even 
Higgs Sa- The Saydjdi coupling reads 



Qal2 = ^(5al(iVl5iV24 + iVi4iV25) + 5,2(iVi5iV23 + iVi3iV25) 



+ Sa3{NnN24 + iVl4iV23)) " V2KSa3N,,N25 
+ ^{92Nll - giNi2){N23Sal - N2iSa2) 

+ \{g2N21 - giN22){N,,Sal - N,,Sa2) (99) 

where Sai denote the mixing components of Sa in the basis (^6|). The Higgs- 
fermion coupling is 

Qaf = i- for an up-type quark, (100) 

V2hi 

Qaf = /- for a down-type quark or a charged lepton (101) 

v2/i2 

where nif is the fermion mass. With the Higgs propagator Da, the invariant 
amplitude reads 

Y, \Ma\' = AQl^2Qlf\Da{s)Wmlo + m^o + m^om^o - s)s, (102) 



spins 
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leading to the following partial width: 



-0 Sa ~0 f n _ ''f^aU^af "V, 






mi -^ xVf) = ',::. ' z:^ia{v, ^a) (103) 



where the phase space integral la depends on 77 and uoa = ms^/m^o: 

4(,...) = /^ ,.^(__,^LH_Z^, (104) 

Finally, we consider the sfermion exchange. Neglecting the fermion Yukawa 
coupling, the x'iffR vertex is 

fi = giV2QfNi2, (105) 

and for the x^f fh vertex one gets 

/, = V2{g2T,fNa - gii^f - Qf)N,2). (106) 

With the sfermion propagator D ?, the invariant amplitude reads [^ 

El-^/P = flfl[Hl-m-ml.)\D-^{t)\' + t^u 

spins 

+2m^om^osDj{t)Dj{u)^ . (107) 

The last term is an interference term. The partial width is then given by 

r(x5^x;//) = ™^//fe^/). (108) 

I?{ri,uj?) is the phase space integral, depending on 77 and u? = m^o/m?: 



Ifijl.Uf) = -r dx- 



'/^'""^/^ 4V {ujjx-iyx 



i-l'7l f^ 2z + ri'' 



+ ^y^ d^J^ ^^ (^yi _ 2.) - 1)^(1 - 2. + 2x) - 1) -^'"'^ 

Now, we turn to the two body decay. The decay rate for the on-shell scalar 
Higgs production of fig. ^ reads 

IGtt 



ml-xISa) = ^^^^^^X{l,v',u;'a)i^ + '2V + V'-U^l)- (110) 



Finally, let us consider the radiative decay X2 ~^ Xil- The analytic approxi- 
mation used in the main part of the paper includes only charged lepton/" right" 
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slepton loops (for details, see |T^,[T^). The different contributions are shown in 
fig. ^. The decay width reads 






X2 



;iii) 



g^^^ is an effective coupling: 

dxxi ~ 



Y^l^^'^f^fhh^^iv^^fj 



f 



;ii2) 



where /i^2 are given in eq. (|105|) and 



IyiV,UJf) 



1 — rj^ Jo 



dx \1 + r] + 



1 — uj'jx 



1 — rjuj'^^x 



log 



(113) 
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Figure Captions 

figure 1: Diagrams contributing to the three body decay X2 ~^ Xiff- 

figure 2: Feynman graph for the two body decay X2 ~^ Xi^a/Pa- 

figure 3: Contributions to the radiative decay X2 ~^ Xil with only charged 
lepton/" right" slepton loops. 

figure 4: Branching ratio of the charged lepton production X2 ~^ Xi^^^^ versus 
the Bino mass m-^o. 

figure 5: Branching ratio of the neutrino production x\ ~^ Xi^^ versus the Bino 
mass 171-^0. 

figure 6: Branching ratio of the jet production X2 ~^ XiQQ {q = ""5 d^ c, s, b) 
versus the Bino mass m^o. 

figure 7: Branching ratio of the real singlet Higgs scalar production x\ ~^ x\Si 
versus log(A). 

figure 8: Branching ratio of the radiative decay x% ~^ Xi7 versus rj. 

figure 9: Total decay width of the bino to singlino transition log(r(x2 ~^ Xi^)) 
versus log(A). 
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